A phase operator for photons 
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We define a unitary phase operator for photons in a single momentum mode. The operator acts 
on a Hilbert space with basis consisting of aU number states in both polarizations. The Susskind 
Glogower operator, E — e'^ , thus defined is unitary and acts as a true ladder operator over all 
the states. It satisfies the commutation relation |^_E, nj = E; consequently, the phase tj) satisfies the 

canonical commutation relation |^n, (^j = i. The eigenstates of the Susskind Glogower operator are 
found. Our model is able to directly account for phase measurements when polarization changes 
cannot be neglected and at the same time gives a well defined phase operator. 

PACS numbers: 42.50.-p, 42.50.St, 03.65.Wj 



I. INTRODUCTION 

The idea of a phase operator for oscillators and fields 
has a long history. Dirac in an attempt to define 
the quantum equivalents of action and angle variables, 
decomposed the annihilation operator for the harmonic 
oscillator as d = EVfi = e^'^\ffi. Here ^ is the phase 
operator and h the number operator of the oscillator. 
The number operator was to be the quantum equivalent 
of the classical action and satisfy canonical commutation 



relations with the phase. 



n, 



As he pointed 



out, this phase operator is not well defined as there is a 
lower bound on the eigenvalues of the number operator. 
In addition, the phase is periodic and defined modulo 27r. 

As almost all measurements on electromagnetic waves 
are measurements of intensity or phase; there has been 
an abiding interest in a phase operator for light too. A 
large number of subsequent suggestions have been made 
for the quantum phase operator. Susskind and Glo- 
gower suggested using the corresponding exponential 
operator as a ladder operator, E\n >= |n — 1 >. As 
the operator E annihilates the vacuum state to give 0; 
it is not a unitary operator. Carruthers and Nieto [1] 
suggested using the sine and cosine operators, which 
could be formed out of the exponential operator and 
its hermitian adjoint but pointed out that they do not 
commute. Carruthers and Nieto also pointed out that 
the quantum number specifying the states was closely 
connected with the spin angular momentum component 
of the beam along the direction of propagation, so 
extending the energy spectrum to negative values would 
solve the unitarity problem. However, they rejected 
this extension as negative energies are not possible. 
This concept of doubling the Hilbert space was revived 
by Newton Q, who also rejected the negative energy 
states but suggested a formal renumbering of the energy 
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state with the states containing n quanta of negative 
energy being labeled — (n + 1). Modifications of this 
procedure were sugg ested by others [5|, 1^, i, [HI. These 
approaches suffer from either lack of unitarity of the 
exponential operator or failure to obtain the canonical 
commutation relations for the number operator with 
the phase operator. An excellent review of the above 
ideas as well as those to be described below was given 
by Lynch 

Pegg and Barnett [1, 0] suggested the most popular 
way round these difficulties. The idea is to truncate the 
oscillator states to a finite number say s. The operators 
in the resulting finite dimensional Hilbert space are well 
behaved, this is extremely useful as a practical tool. 
However, questions have been raised [8l-[l0| of whether 
such a weak convergence to the limiting infinite dimen- 
sional operator gives the same results as the calculations 
directly using infinite dimensional operators. In view of 
these difficulties, Noh, Fougeres and Mandel [Ulil even 
suggested that a unique phase operator should not be 
sought. They suggest that we should only use an opera- 
tional definition of phase based on experiment performed. 

With the recent development of quantum tomographic 
methods to study entanglement and deduce the Wigner 
distribution, the phase or more accurately the phase 
difference has again become very important [Tsj . There 
has been many subsequent efforts to define a phase 
operator both for photons and for the closely related 
problem of the Harmonic oscillator (T6l - l22j . Many of 
these, especially the work by Ban, Luis and Hradil, use 
two mode operators which are especially well suited to 
study parametric processes and higher order interference 
but are also used to cover the case of two orthogonal light 
polarization modes. However as pointed out by Bjork 
et. ah, these too are not well defined in the vacuum 
state and indeed require a separate definition of the two 
mode vacuum state. There is also a related difficulty, 
the two boson representations used are isomorphic to the 
angular momentum group, the Cartesian components 
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of which do not commute. It is known that the the 
Cartesian components of the hght polarization need to 
commute t o g ive a definition consistent with Maxwell's 
equations (23l [23 |. Note that these quanta have been 
measured and shown to add algebraically, with quanta 
of azimuthal angular momentum [2^. 

Sperling and Vogel recently pointed out that re- 
moval of a right circularly polarized photon is the same 
as the addition of a left circularly polarized photon. They 
then suggest treating the left circularly polarized photons 
as antiparticles of the right circularly polarized photons 
and used this to construct a unitary Susskind Glogower 
operator E. The difficulty with this theory is that there is 
a difference in angular momentum between the vacuum 
for left circularly polarized photons and that for right 
circularly polarized photons. 

We suggest a modification of the Sperling- Vogel pro- 
cedure, based on two requirements. First, the quantum 
phase being an angle, the corresponding momentum must 
be an action variable. The number operator acting on a 
Fock state then measures the number of action quanta in 
the state. No experiment directly measures the action, 
but the number of action quanta is equal to the num- 
ber of azimuthal angular momenta in that state, which 
is measurable so we have a suitable proxy. The essence of 
our suggestion is that action and phase form a conjugate 
pair and energy and time form another conjugate pair 
and treating them as freely exchangeable may be part of 
the reason for our difficulties. 

Second, by passage through a phase plate, it is possi- 
ble to transform a state of definite momentum and po- 
larization "a mode" , into one of the orthogonal polariza- 
tion, "an orthogonal mode" , while preserving the phase 
or changing it trivially. So the Hilbert space relevant to 
phase measurements is that spanned by both the polariza- 
tion modes taken together. 

The Hilbert space we consider is a tensor product of 
a two internal dimensional space of polarization and the 
Fock space of a single polarization mode. This agrees 
with the representations found by Wigner, [2^, (23] for 
massless relativistic particles. As Bargmann and Wigner 
showed the appropriate relativistic representation for 
photons is two dimensional, this was elaborated later 
|30| to show that it is indeed the natural representation 
to take into account all polarization phenomena. As the 
detailed calculations of Kim et. al. show the operators 
for the commonly used devices like polarizers, phase 
plates are particularly simple in this representation 
being constant matrices. The action of such operators is 
similar on states with the same polarization but different 
number of photons. Such state independent, passive 
operators are far more difficult to construct in the two 
mode theories of polarization. 

This method permits a natural extension to cover the 
negative integers. If the number of right circular polar- 
ized quanta present is denoted by > then the num- 



ber of left circular quanta present is obviously n_ < or 
vice versa. Newton [^J, first suggested the use of such a 
two dimensional representation, but in terms of energy 
eigenstates, consequently had to suppress the negative 
energy states. Our interpretation though mathematically 
very similar is in terms of polarizations, where the inter- 
nal two space is essential. Thus our model is able to 
directly account for phase measurements when polariza- 
tion cannot be neglected while also giving a well defined 
phase operator. 

In the next two sections, the Hilbert space and the 
phase operator are defined and the operator is shown 
to be Hermitian while the corresponding exponential 
Susskind Glogower operator is unitary. We then find the 
eigen states of the Susskind Glogower operator and show 
they are orthogonal. All the previous results cited above 
can be recovered by restricting our space to states of a 
single polarization. While we outline our results using 
the two opposite circular polarizations, they hold in any 
polarization basis. 



II. THE STATE SPACE 

Consider an electromagnetic field propagating in a vac- 
uum, the states are completely specified by the momen- 
tum and the spin angular momentum component (or po- 
larization) along the direction of propagation [23, |28| . 
The positive frequency components of a single momen- 
tum mode of the radiation field are usually written as. 



E{f,t) oc a+e*^*^-*'-"*) + a^e^^^'^"""*). 



(1) 



Here a_)- is the annihilation operator for right circularly 
polarized light, and a_ is the annihilation operator for 
left circularly polarized light. These are subject to the 
commutation relations. 



a-|_, a_ I = a+, a 



At 



= 1, 
= 0. 



(2) 



There is a modified representation, valid for zero mass 
particles, found by Wigner [l^,!!^], which we give below. 
The two opposite circularly polarized beams are distin- 
guished by their helicity. As there are two states for a 
given momentum k, we write these as a single state with 
an internal degree of freedom. 



\n+ > 
\n- > 



(3) 



They are eigenstates of the Pauli Lyubanski helicity 
operator Pp l , 



PpL 

PpL 







" |n+ > " 





= + 













n_ > 




|n_ > 



(4) 
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III. THE PHASE OPERATOR 

To incorporate the polarization directly into our de- 
scription, we note that the commutation relations for the 
annihilation and creation operators imply that if |n > is 
an eigenstate of the number operator with eigenvalue n 
then 



a^a{d\n >) = (n — l){a\n >), 
d'^ a{a'' \n >) = (n + l){d''\n >). 



(5) 



Both are eigenstates with eigenvalues n—1 and n+1. The 
commutation relations impose no further restrictions on 
n. The series termination requirement yields the fact that 
n must be an integer. However the set of all nonnegative 
integers as well as the set of all non positive integers both 
satisfy equation JS]). 



it 



/nZIn I - 1 >, 



a+\n+ > = A/n+ + 1 |n+ + 1 >, 
a+|0>=0, ^ n+ = 0,+l,+2,+3. . 
a_|n_ > = ^|n_| |n_ - 1 >, 



aL|n_ > = v^l"- + 1| \n- + 1 >, 
1>=0, =^ n_ = -1,-2,-3.... 



(6) 



(7) 



In the subspace n+ > 0, the a+ acts as a annihilation 
operator, while in the subspace n_ < the annihilation 
operator is a[_ . The choice of quantum number sequence 
for the negative quantum numbers is the same as that 
used by Newton, [3,[2^ and serves to emphasize the sim- 
ilarity to the familiar case of positive quantum number 
states. It should be emphasized that the state | — 1 > 
is the vacuum state for the lower polarization. Thus the 
state space is spanned by the basis, 



\n > 




n > 0, and 





In > 



n <0. 



(8) 



Note that though circular polarizations are the natural 
basis of the hclicity operator, the above representation 
can be used with any other polarization basis. 

Any of the operators in a definite polarization sub- 
space can now be extended to the doubled space of both 
polarizations as. 



a+,Pj 



PL 



-0, 



PpL =0. 



For any given n, the state n = n+ and n 



(9) 



I- + 1| 

contain the same number of photons but have opposite 
polarizations. It is important to note that both a+ and 
d}_ are operators which reduce the number of photons in 
their subspace. That is their repeated application will 
cause all states with a finite number of photons to be 
annihilated. 

The method can be used to extend the operation of 
all operators in either subspace to cover the whole of the 



doubled space. So far however, the extension is purely 
formal as we have no operator which connects states la- 
beled by n+ and those labeled by n_. Consider a new 
operator. 



{d\)d, = |0 X 0|, 
dv{dl) = I - 1 >< -1|, 



1 >< 0| 



I + PPL 



(10) 



This operator causes the vacuum of the right circular 
polarized states to become the vacuum of the left circu- 
lar polarized states. A physical realization of this is the 
passage of a polarized vacuum state through a half wave 
plate. The hermitian adjoint couples the subspaces in 
the reverse direction. We use equations ^ and pUj) to 
define a lowering operator which acts on every state in 
the doubled Hilbert space i?. 



d+ 
a„ d- 



(11) 



The commutation relations for the modified annihilation 
and creation operators. 



PpL - 10 >< 0| 

PpL + \-l>< 



-II 



It is instructive to find the matrix elements of equation 
(IT^ between all the states. Obviously, the non-diagonal 
elements all vanish, while the diagonal elements yield 



< n+\ [a™,aJ,J \n+ > = 1, 

< n_| [a„,aj„] |n_ > = -1, 

< -1| [a,„, aL] I - 1 >= < 0| [a™, d^ |0 >= 0. (13) 

We see that the usual commutation relations between 
annihilation and creation operators hold, with the very 
interesting exception that the annihilation and creation 
operators commute with each other between the vacuum 
states |0 > and | — 1 >. For a harmonic oscillator, it is 
the commutation relations eq. ([S]), which give rise to the 
energy difference between the states, so the commutation 
relation (jl3p guarantees that there is no difference in en- 
ergy between the two vacuum states. As the particles are 
bosons, the vacuum state is then 





"|o>- 














+ 


.l-l>. 





(14) 



Following Newton we write the number operator as. 







d^dl-l 



(15) 



This choice is governed by the fact that our operator 
must count the number of photons or the energy. 
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Next we use equation (fTTj) to define a Susskind Glo- 
gower operator, 



Em — 







a„ a- 



(16) 



This operator has all the required properties, 
E„ 



\n+ > 




|n+ - 1 > 




En 





|n_ > 





|n_ - 1 > 



— En 



Along with 



Efn Ejji — Em EfYi 1 • 



(17) 



(18) 



As the Susskind Glogower operator is unitary from 
equation ([18]), we can write Em = e"'^™ and see that 
the hermitian phase operator <f>m satisfies the canonical 
commutation relations 



(19) 



Consequently, we have the canonical representation n„ 



''dip 



for the number operator in a basis of cj) eigenstates. 



IV. THE EIGENSTATES 

We note that the equation (|17p becomes a differen- 
tial equation in the canonical representation. As the 
Susskind Glogower operator is unitary, we can see that 
all eigenstates of this operator must be periodic with pe- 
riod 2tt. Then the Floquct theory of periodic differential 



equations guarantees that every solution of equation p7)) 
must be of the form, e^'^^e*^"'^, where 00 is arbitrary and 
n is an integer. As is well known, these functions form a 
complete basis set for (f) in the interval -I- (/)o to 27r + (/)o. 
So the eigenvalue equation for Em is 



(20) 



with the solution 



>= 



'2ti 



E 

n=0 

E 






g«("+l) 



|n > 
\n > . 



(21) 



It is interesting to note that the two vacuum state con- 
tributions differ by a phase. 



e 2 



|0> 






- 1 > 



(22) 

This suggests that we may be able to distinguish between 
them using the polarization shifts at reflection. Further, 
in agreement with the requirements from energy quanti- 
zation, the action in the vacuum states is also i. 

By restricting our space to states of a single polariza- 
tion, we get back the phase operators defined by other 
authors. We outlined our results using the two opposite 
circular polarizations, but they hold in any polarization 
basis. Detailed calculations including the effects of polar- 
ization separation by Fedosov Imbert or Goos Hanchen 
shifts at reflection from dielectrics will be presented else- 
where. 

CP wishes to acknowledge support from the U.G.C. of 
India. 



[1] P.A.M. Dirac, Proc. Roy. Soc. London Ser. A 114, 243 [12 

(1927). [13 

[2] L. Susskind and J. Glogower, Physics 1, 49 (1964). 

[3] P. Carruthers and M.M. Neito, Rev. Mod. Phys. 40, 441 [14 
(1968). 

[4] R.G. Newton, Ann. Phys. (N.Y.) 14, 327 (1980). [15 
[5] J.C. Garrison and J. Wong., J. Math. Phys. 11 2242 

(1970). [16 

[6] S.M. Barnett and D. T. Pegg, J. Phys. A 19, 3849 (1986). [17 
[7] D.T. Pegg and S.M. Barnett, Europhys. Lett. 6, 483 

(1988). [18 

[8] A. Luis and L.L. Sanchez-Soto, Phys. Rev. A 48, 752 [19 

(1993). [20 
[9] J. Bergou and B-G. Englert, Ann. Phys. (NY) 209, 479 

(1991) . [21 
[10] A. Vourdas, Phys. Rev A45, 1943 (1992). [22 
[11] V.N. Popov and V.S. Varunin, J. Mod. Opt. 39, 2121 

(1992) . [23] 



R. Lynch, Phys. Reps. 256, 357 (1992). 

J.W. Noh, A. Fougeres and L. Mandel, Phys. Rev. A45, 

424 (1992). 

J.W. Noh, A. Fougeres and L. Mandel, Phys. Rev. A46, 
2840 (1992). 

S.L. Braunstein and P. van Loock Rev. Mod. Phys. 77, 
513 (2005). 

Masashi Ban J. Opt. soc. Am. B 9 (1992) 1189. 

M. Daoud and E. H. El-kinani, Phys Letts, A (2011) 

doi:10.1016/j.physleta.2011. 03.018. 

Xian-ting Liang, Phys. Lett. A 273 (2000) 218. 

Z. Hradil, Phys. Rev. A 47 (1993) 93. 

P.L. Garcia de Leon and J. P. Gazeau, Phys. Letts., A 

361 (2007). 

Shunlong Luo, Phys. Letts. A 275 (2000) 165. 
G. Bjork, J. Soderholm, L.L. Sanchez-Soto, A.B. Klimov, 
I. Ghiu, P. Marian and T. A. Marian, arXiv:1 00O935l 
Chun-Fang Li, Xi. Chen Arxiv: 1011.3608 



5 



[24] S.J. van Enk and G. Nicnhuis, J. Mod. Opt. 41(1994) 

963, and EuroPhys. Letts. 25 (1994) 963. 
[25] N.B. Simpson, K. Dholakia, L. Allen and M.J. Padgett, 

Opt. Letts., 22 (1997) 52 
[26] J. Sperling, and W. Vogel, Arxiv:09073349. 
[27] "Coherence and Quantum Optics", L. Mandel and E. 



Wolf, Cambridge University press, Cambridge (UK). 
[28] V. Bargmami and E.P. Wigner, Proc. Nat. Aead. Sei. 
34,211 (1948). 

[29] Y.S. Kim and E.P. Wigner, Phys. Rev. A38, 1159 (1988). 
[30] D. Han, Y.S. Kim and M.E. Noz, Phys. Lett. A219, 26 
(1996). 



